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Abstract 

The multiple- quantum operator algebra formalism has been exploited to 
construct generally an unsorted quantum search algorithm. The exponen- 
tial propagator and its corresponding effective Hamiltonian are constructed 
explicitly that describe in quantum mechanics the time evolution of a multi- 
particle two-state quantum system from the initial state to the output of 
the unsorted quantum search problem. The exponential propagator usually 
may not be compatible with the mathematical structure and principle of the 
search problem and hence is not a real quantum search network, but it can be 
further decomposed into a product of a series of the oracle unitary operations 
such as the selective phase-shift operations and the nonselective unitary op- 
erations which can be expressed further as a sequence of elementary building 
blocks, i.e., the one-qubit quantum gates and the two-qubit diagonal phase 
gates, resulting in that the decomposed propagator is compatible with the 
mathematical structure and principle of the search problem and hence be- 
comes a real quantum search network. The decomposition for the propagator 
can be achieved with the help of the operator algebra structure and symme- 
try of the effective Hamiltonian, and the properties of the multiple-quantum 
operator algebra subspaces, especially the characteristic transformation be- 
havior of the multiple-quantum operators under the z-axis rotations. It has 
been shown that the computational complexity of the search algorithm is 
dependent on that of the numerical multidimensional integration and hence 
it is believed that the search algorithm could solve efficiently the unsorted 
search problem. An NMR device is also proposed to solve efficiently the un- 
sorted search problem in polynomial time. 



1 



1. Introduction 

Quantum computation is a cross discipline among mathematics, quantum 
physics and information science. It obeys both the quantum mechanical laws 
and the mathematical principles [1-5]. It has been discovered that quantum 
computers can solve efficiently certain problems in polynomial time that can 
not be solved efficiently by any classical digital computers [3, 6-11]. The 
famous examples include the prime factorization [10] and the quantum simu- 
lation [3, 11] whose polynomial-time quantum algorithms have been discov- 
ered. However, in practice these problems including the prime factorization 
and the quantum simulation are rather specialized. An important question 
in quantum computation therefore arises whether or not quantum computers 
can solved efficiently a general NP problem in polynomial time [10, 12, 13]. 
It is well known that NP-problems are hard for any classical computation. In 
the past several years this question has been discussed extensively but has re- 
mained largely ignored [12, 13]. On the other hand, it has been found that a 
quite broad class of problems such as search and optimization problems can 
be speeded up quadratically by quantum computation [14]. The unsorted 
search problem is really a hard problem in classical computation. Assume 
that there is a large unsorted database T — (0, 1, ...,N — 1) {N = 2"), in 
which only one of these elements satisfies the function f{s) — 1, but /(r) = 
for any other element r {r ^ s). Now one wants to find the target element s. 
If a classical digital computer is used to search for the target element, one will 
need to examine elements of the database in the worst case and an average 
of A^/2 elements before finding the desired element s. However, Grover [14] 
has showed recently that if a quantum computer is used one needs to exam- 
ine only ^/N elements around to find the target element. It has been proven 
that the Grover algorithm is the optimal quantum search algorithm so far 
[15]. Grover and his coworkers has shown further how his search algorithm 
can speed up quadratically almost any other quantum algorithms [16, 17]. 
However, the Grover algorithm is not an efficient quantum search algorithm 
and therefore can not solve efficiently a general NP-problem in polynomial 
time [12]. 

One of the most important characteristic features for quantum compu- 
tation different from classical counterpart is that quantum computers can 
have the ability of the massive parallel processing in computation [4]. In 
principle, a quantum computer can offer the possibility for solving efficiently 
a general NP problem in polynomial time by virtue of the massive quantum 
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parallelism. However, the quantum computational output usually can not 
be obtained correctly and directly in polynomial time due to the limits of 
quantum mechanical measurement. This is the reason why there exists the 
question whether a quantum computer can solve efficiently a general NP- 
problem or not. As a consequence of the Grover quantum search algorithm 
[14] , it has been shown that the correct quantum computational output may 
be obtained ahci\^ iterations and therefore, a general NP problem may 
be speeded up quadratically by quantum computation [12]. It has been be- 
lieved extensively that quantum computers can solve efficiently a general 
NP-problem by virtue of the massive quantum parallelism, although so far 
only few polynomial-time quantum algorithms have been discovered to solve 
efficiently some special hard problems [4, 10, 13]. 

The multiple-quantum operator algebra formalism has been proposed to 
describe quantum computational process [18, 19]. It has been exploited ex- 
tensively to design the quantum computational network of a known quantum 
algorithm [19], construct elementary building blocks of quantum computa- 
tion [19, 20a], and prepare the effective pure states in the NMR quantum 
computation [20b]. In this paper the multiple-quantum operator algebra for- 
malism has been used to construct an unsorted quantum search algorithm 
and its quantum computational networks that obey the quantum mechani- 
cal laws and are compatible with the mathematical structure and principles 
of the unsorted search problem. This quantum search algorithm is different 
from the Grover's one in that both two algorithms have different propagators. 
Its computational complexity is dependent on that of the numerical multi- 
dimensional integration. Therefore, it is believed that the quantum search 
algorithm could solve efficiently the unsorted search problem on an oracle 
universal quantum computer. This quantum search algorithm is constructed 
with two families of elementary unitary operations, that is, the nonselec- 
tive unitary operations and the oracle unitary operations, e.g., the selective 
phase-shift operations, all these operations can be further decomposed into 
a product of a series of elementary building blocks such as one-qubit quan- 
tum gates and the two-qubit diagonal phase gates [19]. Here assume that 
the quantum search algorithm runs on a quantum computer with a quantum 
system consisting of n two-state particles such as n nuclear spins with the 
angular momentum quantum number 1=1/2 and one can manipulate at will 
each individual two-state particle of the quantum system by an external field 
such as an electromagnetic field. Also assume that any decoherence effects 
are ignored in the quantum system. 
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2. The effective Hamiltonian of a quantum system and the con- 
struction of quantum algorithms 

Benioff was the first time to use quantum mechanics to describe the re- 
versible computation process on classical Turing machines by understanding 
the corresponding relationship between the reversible computation process 
and the fact that time evolution of an isolated quantum system is reversible 
dynamic [1, 2]. He showed that quantum mechanical models were computa- 
tionally as powerful as the classical Turing machines. Feynman was the first 
person to conjecture that quantum mechanical models might be more pow- 
erful than any classical computers in simulating quantum processes [3]. His 
universal quantum simulator could efficiently simulate any quantum dynam- 
ics of a quantum system whose Hamiltonian consists of any local interactions 
[11]. Deutsch formalized the concept of the universal quantum computer and 
has showed that quantum Turing machines could be more powerful than the 
classical counterpart from a computational complexity point of view [4] . He 
has also developed the quantum circuit model of quantum computation [5]. 
Yao has showed further that the two models of quantum computation, i.e., 
the quantum Turing machine and the quantum circuit model, are polynomi- 
ally equivalent to each other [21]. Therefore, there are two ways of thinking 
about quantum computation [4, 5, 8]. One way is to think of it as the re- 
versible computation on a quantum Turing machine, and another is that 
quantum computation can be thought of as the unitary time evolution of a 
quantum system from the input state to the output. The unitary time evo- 
lution can be described by a unitary transformation that acts on the input 
quantum state in a quantum system. In quantum mechanics, time evolution 
of a quantum system may obey merely quantum mechanical laws, e.g., the 
Schrodinger equation. Time-evolutional propagator of a quantum system 
characterizes completely the unitary dynamical behaviors of the quantum 
system. However, in quantum computation time evolution of a quantum sys- 
tem during quantum computing from the input quantum state to the output 
is also subjected to the mathematical structure and principle of the prob- 
lem solved by a given quantum algorithm running on the quantum system 
in addition to these quantum mechanical laws [19]. Therefore, the form of 
the propagator is constrained by the mathematical structure and principle 
of the quantum algorithm. This propagator may be considered generally 
as a unitary transformation that may be taken as an exponential unitary 
operator U{t) — exp{—iHt), where the operator H is the effective Hamil- 
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tonian of the quantum system subjected to the quantum algorithm. This 
shows that the form of the effective Hamiltonian H is also constrained by 
the mathematical structure and principle of the quantum algorithm. Then 
the effective Hamiltonian H could also really characterize the mathemati- 
cal structure and properties of the quantum algorithm such as the quantum 
computational complexity. For example, provided that the effective Hamilto- 
nian H consists of any local interactions in the quantum system, there could 
be a quantum network that can simulate efficiently time evolution of the 
quantum system subjected to the quantum algorithm [3, 11]. Obviously, the 
quantum network is always compatible with the quantum mechanical laws. 
Is the quantum network also compatible with the mathematical structure 
and principle of the problem solved by the quantum algorithm? Evidently, 
if the quantum network is designed according to the mathematical struc- 
ture and principle of the problem it should be a real quantum computational 
network to solve the problem, otherwise it can not be thought of as a real 
quantum network to solve the problem. Now, if the problem solved by the 
quantum algorithm is an NP-problem this real quantum computational net- 
work can solve efficiently the NP-problem because it can simulate efficiently 
time evolution of the quantum system subjected to the quantum algorithm. 

A quantum computation is a unitary time-evolutional dynamical process 
subjected to a given quantum algorithm from the input quantum states to 
the output in a quantum system. Now one is given a mathematical problem 
that needs to be solved on a quantum computer. Assume that there are a 
number of quantum algorithms to solve the same problem. In practice, one 
may also design a number of quantum computational networks for a given 
quantum algorithm to solve the same problem. Consider the special case 
that the initial state and the output in the quantum system are fixed for 
a given mathematical problem to be solved on the quantum computer. For 
example, for the unsorted quantum search problem the initial state is usually 
considered as the superposition and the output state is the target state in a 
quantum system. It is well known that there are a number of unitary time- 
evolutional pathways from the fixed input state to the output in a quantum 
system. Each such pathway is governed by the quantum mechanical laws 
and described by a time-evolutional propagator or its corresponding effective 
Hamiltonian. There may be a unitary dynamical process of the quantum 
system subjected to the quantum algorithm to solve the problem among all 
these unitary time-evolutional pathways, and this process is characterized 
completely by the effective Hamiltonian during running the quantum algo- 
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rithm. This unitary time-evolutional process obeys not only the quantum 
mechanical laws but also the mathematical structure and principle of the 
problem. If there are a number of quantum algorithms to solve the same 
problem with the fixed input state and the output, there is also a unitary 
dynamical process corresponding to each such quantum algorithm that obeys 
both the quantum mechanical laws and the mathematical principles of the 
problem. Now, assume that the problem is an NP-problem in classical com- 
putation, is there an efficient quantum algorithm to solve the problem among 
these quantum algorithms? If the efficient quantum algorithm exists how to 
find it and construct its quantum computational network? These problems 
have not be solved and reminded largely ignored to date. One the other 
hand, there already exists a simple scheme to find a quantum algorithm to 
solve a mathematical problem when there is a classical algorithm to solve 
the same problem. This simple scheme is that a quantum algorithm can be 
obtained from its corresponding classical algorithm. The computational net- 
work of the classical algorithm usually consists of a sequence of irreversible 
classical logical gates. It can be translated into a quantum algorithm with 
at least the same computational power by simply replacing the irreversible 
classical logic gates with the corresponding reversible quantum gates accord- 
ing to the Bennett's suggestion [22]. However, this simple scheme usually 
may not be available for finding an efficient quantum algorithm to solve a 
general NP-problem. In this paper the multiple-quantum operator algebra 
formalism is exploited to design a new quantum algorithm and construct its 
quantum computational network for a given mathematical problem such as 
the unsorted search problem [18], where it is assumed that the input state 
and the output of the quantum system for the problem to be solved are given 
in advance. Then the total time-evolutional propagator that transforms uni- 
tarily the input state to the output can be constructed explicitly. In general, 
the propagator obeys the quantum mechanical laws but usually is not com- 
patible with the mathematical structure and principle of the problem. Then 
it is not a real quantum algorithm to solve the problem. However, the prop- 
agator may be further decomposed into a sequence of the quantum circuit 
units which are compatible with the mathematical structure and principle of 
the problem with the help of the properties of the multiple-quantum opera- 
tor algebra spaces [18, 19] and the operator algebra structure and symmetry 
of the effective Hamiltonian. Finally these quantum circuit units are fur- 
ther decomposed into a product of a series of elementary building blocks, i.e, 
the one-qubit quantum gates and the two-qubit diagonal phase gates [19]. 
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Such constructed quantum network is obviously governed by the quantum 
mechanical laws and compatible with the mathematical structure and prin- 
ciple of the problem to be solved. Therefore, it becomes a real quantum 
algorithm to solve the problem. Obviously, once the propagator is decom- 
posed into a product of a polynomial number of elementary building blocks 
the constructed quantum network is really a polynomial-time quantum com- 
putational network to solve efficiently the problem no matter whether the 
problem is an NP-problem or a polynomial-time problem in classical compu- 
tation. As an example, a new quantum search algorithm to find the marked 
element in an unsorted database is explicitly constructed with the help of the 
multiple-quantum operator algebra formalism. Its quantum computational 
network consists of the two types of elementary quantum circuit units, that 
is, the oracle quantum unitary operations, e.g., the selective phase-shift oper- 
ations, and the nonselective unitary operations, i.e., the oracle- independent 
quantum unitary operations. 

3. The nonselective unitary operations 

A nonselective unitary operation is an oracle-independent unitary opera- 
tion that acts on every two-state particle of a quantum system symmetrically, 
that is, all the two-state particles in the quantum system are indistinguish- 
able and symmetrical with respect to the unitary operation. The nonselective 
unitary operations are independent of the marked state to be searched for 
in the quantum system and can be implemented on a quantum computer 
without knowing in advance any state of the quantum system. They are 
also independent of whether the database under search is sorted or unsorted. 
Therefore, this type of unitary operations are compatible with the mathe- 
matical structure and principles of the search problem and can be used to 
build quantum computational networks of a quantum search algorithm. The 
Walsh-Hadamard transformation W and the diffusion transform D in the 
Grover algorithm [14] are typical nonselective unitary operations. Besides 
the two unitary operations there are a number of other nonselective unitary 
operations. As an example, two types of general nonselective unitary opera- 
tions are given below, which may be encountered in present quantum search 
algorithm. One type of the nonselective rotation operations that are applied 
to all the two-state particles in a quantum system symmetrically are defined 
by 

Rp{e,m)^eK^{-ieF^) {m ^ 1,2, ...■,p = x,y, z). (1) 
These nonselective unitary operations are constructed with the mth power of 
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the symmetrical Hermitian operator Fp defined by (in a spin- 1/2 language) 

n 

Fp=Y. hp (2) 

k=l 

where the magnetization operator I^p — \(Tkp {p = x, y, z), is the Pauli s 
operator of the kih. two-state particle of the quantum system. Actually, the 
Walsh-Hadamard transform W can be decomposed into a product of the 
nonselective unitary operations Rp{6, 1) [19]: 

W = exp(m7r/2) exp(-i7rF^) exp(-i|Fy). (3) 
Another type of nonselective unitary operations are defined by 

n 

T{9, a, f3,, f3y, 13,) = exp(-^^ A^) (4) 

fe=i 

where the Hermitian operator Aj. of the A;th two-state particle of the quantum 
system may be generally chosen as (in a spin- 1/2 language) 

Ak = aEk + 213 Jkx + 'ipyhy + 213 Jkz (5) 
and the parameter vector {j3p\ is a vector with unit magnitude. The real 
parameters a and 13 p (p = x,y,z) are independent of any index k, indi- 
cating that the unitary operator T{9,a, (3^., Py, (3,) is a nonselective unitary 
operation. Actually, the nonselective unitary operation can be converted uni- 
tarily into simple nonselective phase-shift operations. One of the nonselective 
phase-shift operations is given by [23, 24] 
Co{/3)=Dtag[e-'^^l,...^.,l] 

- exp[-i/?(iEi + lu) ®{\E2 + hz) <S> - (8)(|^n + hz)]- (6) 
When a — 1 the nonselective unitary operation T{9,a, (3^., I3y, (3,) (4) can 
be transferred into the nonselective phase-shift operation Co(/3). Another is 
defined by 

S{I3) = exp(-2/5 X 2Ji, (8) 2/2, (8) ... (2) 24,). (7) 
When a = the unitary operation T{9,a, P^, Py, P,) (4) can be converted 
unitarily into the nonselective phase-shift operation S{P). In particular, the 
general diffusion transform can be constructed with the nonselective unitary 
operations mentioned above: 

D{9) = -WCoi9)W. (8) 
It can prove easily that D{9) = —E + (1 — e^"')P, where E is unity operator 
and the project operator Pij — 1/N, for all i, j. Clearly, the diffusion trans- 
form in the Grover algorithm D — D{tt) [14]. 

4. The selective unitary transformation and the oracle quantum 
unitary operation 

The selective unitary operations are related only to the marked state |s) 
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that is to be searched for in the quantum system. A type of particularly im- 
portant selective unitary operations in an unsorted quantum search problem 
are the selective phase-shift unitary operations. For example, the selective 
phase inversion operation Cs for the marked state in the Grover search algo- 
rithm [14] is a typical selective unitary operation. This type of the selective 
phase-shift operations can be defined as the diagonal unitary operator Cs{9) 
(s 7^ 0, — 1) that has diagonal unitary representation in usual quantum 
computational basis: 

Cs{e) = Diag{l, 1, e-'\ 1, 1} = exp[-^^Ej (9) 
where [Cs(^)]ss — e~^^ only for the diagonal index s, and unit for any other 
diagonal index t ^ s. In particular, the selective phase inversion operation 
Cs = Cs{7i) [14]. As an exception, Co{9) and C]\f^i{9) are two nonselective 
phase-shift operations. The diagonal operator Ess of Eq. (9) can be expressed 
as 

Ess^Diag{0,...,0,l,0,...,0} 

= (|£;i + af7u)0(i£;2 + a|72.)(2)...0(|£^n + <4z) (10) 
where is the 2 x 2-dimensional unity operator and a| = ±1, A; = 1, 2, n. 
It is easy to see that 

'^E, + alh,^(^l ° j,ifa| = +l;i£;, + a^7,,= (^° 5j,ifa| = -l. 

When the selective phase-shift operation Cs{0) acts on an arbitrary computa- 
tional basis |r) a phase shift of exp(— i^) is generated if and only if |r) = |s), 

Cs{9)\r)=exp{-i9Srs)\r). (11) 
Obviously, given an n-dimensional unity-number vector {af , ag, a* } of the 
mared state |s) the diagonal operator Ess is determined uniquely through 
Eq.(lO) and vice versa. Therefore, the unity-number vector {af , 03, a^} 
characterizes completely the diagonal operator Ess ^-nd also the selective 
phase-shift operation Cs{9). 

In general, a selective phase-shift unitary operation may be defined as a 
diagonal unitary exponential operator 

Gs{9) = exp(-ie^,) 
where the diagonal Hamiltonian Hg = ^^^(af , .... ofj is dependent only 
upon the unity- number vector {a'l} of the marked state \s). The diago- 
nal Hamiltonian Hs belongs the /ongitudinal magnetization and spin order 
(LOMSO) operator subspace of the Liouville operator space of the two- 
state quantum system and can be generally expressed as a sum of the base 
operators of the subspace [18] 
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Hs — 0,Q+ J2 ^k^kz+ J2 Jkflhzllz+ Jklm^hzllzlmz + ■■■ (12) 

fc=l k<l k<l<m 

where all the coefficients {^2^, J^;, J^;^, ...} are dependent only on the unity- 
number vector {of, 03, a^}. 

A selective unitary operation is a black-box operation in a quantum search 
problem because this operation is dependent on the marked state whereas 
the marked state needs to be searched for in the quantum system. There- 
fore, the selective unitary operations could be implemented only on an oracle 
quantum computer. A quantum oracle may be defined as a device that, 
when called, applies a fixed unitary transformation Uo to the current quan- 
tum state |r) of the quantum system, replacing it by Uo\r) [12]. There are 
some requirements on the unitary transformation Uo in the quantum search 
problem. The unitary transformation Uo is a selective unitary operation or 
can be expressed as a sequence of the selective unitary operations and non- 
selective unitary operations. The effective Hamiltonian of a quantum oracle 
corresponding to the oracle unitary operation Uo{9) — eyi'p{—i6Ho) is there- 
fore expressed as the form 

Ho = Ho{al,a^,...,a'J. (13) 
Another natural restriction to impose upon Ug may be that Ug is periodic, 
that is, C/^ = r is a known integer, so that the effect of an oracle call 
can be undone by further r — 1 calls, two oracle calls undone by further 
r — 2 calls, and so forth, on the same oracle. The key property of a quan- 
tum oracle is its block-box nature [7]. Therefore, an oracle unitary operation 
could be implemented only on an oracle quantum computer. There are a 
variety of oracle unitary operations in a quantum search problem. For ex- 
ample, the selective phase inversion operation Cs{7i) is chosen as an oracle 
unitary operation in the Grover algorithm whose effective Hamiltonian Egg 
is given in Eq.(lO) and the phase angle 9 = ir [14]. As a generalization of 
the selective phase inversion operation Cs{n) one can choose conveniently the 
selective phase-shift operation Cs{0) as the oracle unitary operation. The or- 
acle unitary operation Cs{9) can be really implemented directly on an oracle 
universal quantum computer and by using this oracle unitary operation one 
can construct the unsorted quantum search network, as shown below. Be- 
sides the selective phase-shift operation Cs{0) the selective unitary operation 
Uop{9) — ex.p{—i9Hop) whose effective Hamiltonian is given by 

Hop = E alhp = {alhp) (g) -^2 ... 

k=l 

+Ei ^{alhp) (S)E3 0...0E„ 
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+... + ^i0...0^„_i0(a^7„p) {p^x,y,z) (14) 
also may be a suitable oracle unitary operation used to construct the quantum 
search network because it may be implemented directly on an NMR quantum 
computer, as shown in Appendix A and C. The effective Hamiltonian Hop is 
traceless and the selective unitary operation Uop{9) of the Hamiltonian Hgp 
can be expressed as 

n 

Uop{e)=U eM-iOaihpl- (15) 

fc=i 

Note that the three oracle unitary operations with different p — x,y,z are 
equivalent to each other by a simple nonselective unitary transformation, for 
example, Uoy{0) = exp(i|Fj,)f/o^(^) exp(-i|Fa;) and Uox{d) = cxp{-i^Fy) 
xUoz{0) exp{i^Fy). Obviously, [Uop{9)Y = E when r9 = Att. The oracle uni- 
tary operation is really single-qubit p-axis pulse applied to all the two-state 
particles in the quantum system and the phase of the pulse applied to the 
kth two-state particle is taken as {—9al). It is well known that single-qubit 
quantum operations are always implemented easily in quantum computation. 
Therefore, this oracle unitary operation is very simple. By exploiting the or- 
acle unitary operation Uop{0) one can decomposed the selective phase-shift 
operation Cs{0) into a product of a polynomial number of the oracle unitary 
operations Uop{0) and the nonselective unitary operations 

Csi9) = cxp{-t^Fy)Uoyi-^)Coie)Uoyi^) expit^Fy). (16) 
Obviously, if the oracle unitary operation UopiO) could be implemented ef- 
ficiently on an oracle quantum computer the selective phase-shift operation 
Cs{0) could be performed in a polynomial time on the same oracle quantum 
computer. There is a question whether the oracle unitary operation UopiO) 
can be expressed as a sequence of the selective phase-shift operations Cs{0) 
and the nonselective unitary operations. This question will be discussed in 
the following section 7. 

How to implement the selective phase-shift unitary operation Cs{0)7 As- 
sume that the quantum system is at an arbitrary state including the marked 
state \s): 

m^\r,S)=E' a.k)[^(|0)-|l))] 

where the ancillary qubit S is at the superposition :^(|0) — |1)). The eval- 
uation of the fTinction f{x) then can be achieved by performing the oracle 
unitary operation Uf on the state |^) [12] 

Uf : I*) ^ ax\x)[^(\0®f(x)) - |l©/(x)))] 
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N-l , , 

= E (-lK(^W|a:)[^(|0)-|l))] 

= a.|^)[^(|0)-|l))]-a,|.)[^(|0)-|l))] 

where the function /(s) = 1 and /(r) = 0, r 7^ s . Obviously, only the target 
state as|s)[;^(|0) — |1))] is inverted and any other state keeps unchanged 
when performing once evaluation of the function f[x). Therefore, perform- 
ing once evaluation of the function f{x) is actually equivalent to applying 
the selective phase-inversion operation Cs(7r) to the quantum system. The 
general selective phase-shift operation Cs{&) can be achieved by the following 
oracle unitary operations 

AT-l AT-l 

[//: E a,|x)|0)|l)^ E a,|x)|0©/(x))|l) 

a;=0 x=Q 

N-l 

V{e): -> E cxpH^5(/(a;),l)]a,|a;)|0e/(a;))|l) 

x-=0 

N-l 

Uf: ^ E a.|x)|0)|l)+exp(-i^)a,|s)|0)|l) 

where the function 5{f{x), 1) = 1 if f{x) — 1; otherwise 5{f{x), 1) = 0, and 
two ancillary qubits arc used in the implementation of the selective phase- 
shift operation Cs{0). Therefore, the oracle unitary operation Cs{0) can be 
expressed as 

Cs{e) = Uf'V{9)Uf (17) 
where V{9) is a conditional phase-shift operation applying only to the two 
ancillary qubits [10, 25]. Note that Uf = 1 the selective phase-shift operation 
can be also expressed as Cs{9) — UfV{9)Uf. 

5. Construction of quantum search networks 

Assume that each usual quantum computational base |r) of the quantum 
system corresponds one-to-one to an element of the search database, and in 
particular, the target element is represented by the marked base \s). A usual 
quantum computational basis can be taken as a Kroncckcr product of the 
common eigenbase of the single-particle spin angular momentum operators 
/| and Ikz in a two-state multiparticle quantum system, for example, 

|s) = \af3...a) — \a) \f3) ... \a). 
In the binary representation the eigenbase \a), \f3) are denoted briefly as |0), 
|1), respectively. In the spinor or vector representation the eigenbase can be 
expressed as \a) = and = In the unity-number representation 
{aj, a^} an arbitrary usual computational basis \s) can be expressed 
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explicitly as 

\s) = (|Ti + a^5i)0(iT2 + a^52)(8)...0(iT„ + <5„) (18) 
where = \a)k + \P)k and 25*^ = \a)k — \P)k- This shows that the unity- 
number vector {af, a*} really characterizes completely the usual quan- 
tum computational basis \s) and vice versa. Obviously, ^T^ + alSk = \a) if 
a| = l; ir, + a|5fe = |/3) if a| = -1. 

A general quantum search algorithm to find the marked state in the quan- 
tum system may be thought of as a unitary time-evolutional process of the 
quantum system from the initial superposition |^') into the marked state \s) 

1 7V-1 

Us: m = ^ E \k) \s) (19) 
ViV k=0 

where the unitary operator Us transforms the initial state to the output of 
the quantum system and may be taken as the quantum computational net- 
work of a quantum search algorithm. There are a number of unitary trans- 
formations Us and different unitary transformations Us may correspond to 
different quantum search algorithms. In the Grover search algorithm [14] 
the unitary transformation Us is taken as a sequence of the 0{-\fN) number 
of the simple unitary transformations: —WCq{t{)WCs{ti)^ where the oracle 
unitary operation is the selective phase- inversion operation Cs(7r). Therefore, 
the Grover search algorithm is a quadratically speed-up unsorted quantum 
search algorithm. Now, a simple unitary transformation Us is constructed 
explicitly that is different from the Grovcr's one and its corresponding effec- 
tive Hamiltonian is local. First of all, there is a unitary transformation that 
transforms unitarily the computational basis |r) to the marked state |s) 

Urs\r) = \s). (20) 
This unitary transformation Urs may be simply constructed by 

Urs=E- Err " E^s + 2/^^ (21) 

where E is the unity operator, the operator Ers is defined by 

and the single-transition operators are defined by [23] 

r^ = l(E +E ) F^^^(E -E ) r'' = ^(E -E \ 
It can prove that the unitary operator Uj-s can be further expressed as the 
exponential form [18] 

t/,, = C,(7r)exp(^7r/;^). (23) 
Therefore, one of the unitary transformations Us of Eq.(19) may be expressed 
as 

Us = UosW = C,(7r) exp(m/o^)I^ (24) 
where the following transformation involved in the Walsh-Hadamard opera- 
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tion W has be introduced [26] 

1 N-l 

E \k) = |0). (25) 

VA/ k=0 

The exponential unitary operator exp(i7r/y*) could be prepared and per- 
formed directly on a quantum computer only when the pair quantum com- 
putational base |r) and \s) are known in advance. However, the basis state 
|s) is the marked state that needs to be searched for in a quantum system. 
Obviously, this operation is not compatible with the mathematical structure 
and principle of the search problem. Then the unitary transformation Us 
of Eq.(24) is not a real quantum search network. In order to construct a 
real quantum search algorithm one needs to find another unitary transfor- 
mation that is entirely equivalent to the unitary operation exp(i7r/°*). This 
unitary transformation is required to be compatible with the mathematical 
structure and principle of the search problem. Obviously, if the unitary op- 
eration exp(z7r/y*) can be expressed as a sequence of the nonselective unitary 
operations and the selective unitary transformations, i.e., the oracle quan- 
tum unitary operations, the unitary transformation Us of Eq.(24) will be 
compatible with the mathematical structure and principle of the quantum 
search problem and hence becomes a real quantum search network. Obvi- 
ously, this quantum search algorithm is different from the Grover's one that 
has a different propagator transforming the initial state to the final state 
of the search problem [17]. In the following it is shown that it is possible 
to express the unitary operator exp(i7r/°^) as a sequence of the nonselective 
unitary transformations and the selective unitary operations, and each of 
these unitary operations can be further expressed as a sequence of one-qubit 
quantum gates and the two-qubit diagonal phase gates. 

Now the single-transition multiple-quantum operator 7^* is expressed as 

^ Yi^ErrQEss — EgsQErr) (26) 

where the operator Q is defined by the matrix elements Qij = 1 for all indexes 
i, j, and can be written in the unity- number representation as the form [20b] 
Q = 2«(iEi + /i,) (1^2 + /2.) ... 0(i^„ + In.)- (27) 
By exploiting the operators E^r, Egg, and Q in the unity- number representa- 
tion one obtains that 

n 

ErrQEgs^® h{alal), (28a) 
fe=i 

EggQE„=® gk{al,al) {28b) 

k=l 

where 
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fkial al) = + alai)E, + l{al + a|)4. 

+ 1(1 - alaDhx + ^i{a'k - ak)hy], 
9k{al al) = + alal)Ek + + a^hz 

+ 1(1 - a^a|)/te - liial - aDhy]. 
In particular, aJJ. = +1, A; = 1, 2, n for the case of \r) = |0) and 

AK,%)-{^+ ifa| = -l' 

^ _ r l^k + hz, if flfc = 1 

^.K,%)-{^-^ ifa| = -l 

where = Ikx ^ihy Then, by inserting Eq.(28) into Eq.(26) the multiple- 
quantum operator can be rewritten generally as 

~ ^mi^m2---^rrn:)i2-^rak+i ^ '^mk+iz)---i2-^mn ^ ^rn„z) (29) 

where the Kronecker product symbol and index order are omitted without 
confusion and here assume that a^. = —1, i = 1,2, ...,k; a^. = 1, i — 
k + 1, ...,n; rrii = 1, 2, n. For example, 

^(/+73+ - /r/3")(|^2 + /2z)... = ®{lE2 + /2.) (8) /s"" - 

-J./r0(i^2+/2.)0/3"0-- 

Equation (29) shows that the operator I^"^ is a k-order multiple-quantum co- 
herence operator (0 < A; < n) [23]. Now one makes a unitary transformation 
on the operator ly'' to convert it into the diagonal operator [20b], 

UI^'U+^lEoo-^E,, (30) 
where the unitary transformation U turns out to be taken as the form 

U exp( i^x2 Im\x---Irni^x-^rn]^^-iz---Imnz) ■ (^l) 

This unitary operator can be further expressed as a sequence of nonselective 
unitary operations and the selective unitary operation Uoy{^\) : 
U = exp(-if Fj,)[/o,(-f ) exp(-if x T^hzhz-InzWoyil) exp(if F^). (32) 
It follows from Eq.(30) that the unitary transformation exp(i7r/°*) can be 
expressed explicitly as 

cxp(i7r/o^) = f/+Co(-7r/2)C,(7r/2)?7 (33) 
where the unitary operation Co(7r/2) is a nonselective unitary operation and 
Cs{9) is the selective phase-shift operation. The unitary operator U (and 
U^) consists of the selective unitary operation [/oi/(±f ) and the nonselective 
unitary operations exp(±i^Fj^) and exp(±i| x 2"'Iizl2z---Inz), as shown in 
Eq.(32). Inserting Eq.(33) into Eq.(24) one writes the unitary transforma- 
tion Us SiS the form 

Us = Cs{7r)U+Co{-7r/2)Cs{7r/2)UW. (34) 
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Now, is this unitary operator Us a real quantum search network? Because 
the selective phase-shift operation Cs{9) can be expressed as a simple se- 
quence of the oracle unitary operation Uoyi;^^) and the nonselective unitary 
operations exp(±i|Fj^) and Cq{6), as shown in Eq.(16), the unitary opera- 
tor Us consists of the oracle unitary operation Uoyi^j) in addition to those 
nonselective unitary operations. Then the quantum network Us oi Eq.(34) 
will be a real quantum search network only when the oracle unitary opera- 
tion f/oy(±^) can be directly implemented on an oracle quantum computer. 
I will show in Appendix A and C how an NMR device can be exploited 
to implement directly and efficiently the oracle unitary operation Uoyi^j), 
then the quantum network Us of Eq.(34) becomes a real efficient quantum 
search network on the NMR device. On the other hand, it has been shown 
in Eq.(17) in section 4 that the selective phase-shift operation Cs{0) can 
be implemented directly on an oracle universal quantum computer [12, 14]. 
Then, the quantum network Us oi Eq.(34) will be a real quantum search 
network when the oracle unitary operation C/o2/(±f ) is expressed explicitly as 
a sequence of the selective phase-shift operations Cs{0) and the nonselective 
unitary operations. Obviously, this real quantum search network Us (34) 
is independent of any specific quantum computer such as an NMR quan- 
tum computer and consists of the selective phase-shift operations Cs{0) and 
the nonselective unitary operations such as the Walsh-Hadamard transform 
exp(iti^ X 2^Iizl2z---Inz)i etc., which can be decomposed further into a 
product of a polynomial number 0{n) of one-qubit quantum gates and the 
two-qubit diagonal phase gates [19], as can be seen in next sections. 

6. The parallel quantum search networks 

As shown in next sections, there is a complex expression for the oracle 
unitary operation Uoy{0) as the selective phase-shift operations Cs{9) and 
the nonselective unitary operations. There are four oracle unitary operations 
Uoy{0) in the quantum network Us of Eq.(34). Therefore, the present quan- 
tum network is quite complicated and is not highly efficient. To simplify it 
one first expands the unitary operator of Eq.(32) as 

U = -j=[E-texp{-t^Fy)Uoy{-^){-ir2''IiJ2z...Inz] (35) 
where the following operator identity has been introduced 

exp(±i7rFp) = (±i)"2"/ip/2p.../np {p = x, y, z). (36) 
With the help of Eq.(35) the quantum network Us of Eq.(34) can be expanded 

as 

Us = \{UaW + U+W - {-iY+^Cs{t:)U,C^{-1)W 



16 



-ar+'Coi-^ww} (37) 

where 

Ua = Co( — f )Cs( — |), 

Ub = expi-tTiF,) exp(zf F,)t/o,(f )C,(f ). 
Equation (37) shows that the quantum network Us is the sum of the four 
unitary transformations: UaW, U+W, Cs{7r)UbCo{-^)W, Co{-^)U^W, indi- 
cating that the quantum network Us could be achieved if one performs in 
parallehsm the four unitary transformations on the same initial superposi- 
tion, respectively, and then sums up coherently their outputs according to 
Eq.(37). Obviously, the parallel quantum network of Eq.(37) is more effi- 
cient with respect with the one in Eq.(34) because one needs to perform in 
parallelism only once the oracle unitary operation Uoy{^). 

7. The transformation between oracle unitary operations 

In a quantum search of an unsorted database the selective phase-shift 
operations Cs{0) can be generally implemented directly on an oracle universal 
quantum computer, as shown in section 4. Therefore, one needs to express 
any other oracle unitary operations such as the oracle unitary operation 
Uop{0) ip = x,y,z) in the quantum network Us of Eq.(34) as a product of 
a series of the selective phase-shift operations Cs{0) and the nonselective 
unitary operations. Obviously, the quantum network Us of Eq.(34) is also 
real quantum search algorithm when the oracle unitary operation UopiO) is 
expressed as a sequence of the selective phase-shift operations Cs{0) and 
the nonselective unitary operations. How to implement directly the oracle 
unitary operation Uop{9) through an NMR device in polynomial time is given 
in Appendix A and C. It is discussed below how to construct the oracle 
unitary operation Uop{9) with the selective phase-shift operations Cs{9) and 
the nonselective unitary operations. 

Now the diagonal operator Egs of Eq.(lO) is expanded as 

n n 

2^E,s = E+ E {2allk,)+ E {2alh,){2atli,) 

k=l l=k<l 

n 

+ E {2ailk,){2atli,){2allm,) + .... (38) 

l=k<l<m 

By using the following spin-echo sequence [23] the even-body interactions 
such as {2allkz){2a\liz), (2aiIkz){2alIi^){2a'^Imz){2a'Jnz), etc., on the right- 
hand of Eq.(38) are cancelled. One therefore has 

Uss = exp{-i6Hss) 

— exp{—i62"'Ess) exp{—i7rFy) exp{i62"-Ess) exp^inFy) 
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= exp[-z2e{E (2a|7fe,)+ E (247fe,)(2af7i,)(2a^7„,) +...}] (39) 

k=l l=k<l<m 

where the effective Hamiltonian 77ss is written as 
OHss = NOEss - ex.p{-i7rFy)NeEssex.p(i7rFy) 

n n 

= 2^{E (2a^7,,)+ E (2a^7fc,)(2af7,,)(2a^7^,) + ...} (40) 

fe=l l=k<l<m 

where N — 2^. Since one can manipulate at will any individual two-state 
particle of the quantum system all the interactions involving only the jth 
two-state particle (j = 1,2, ...,n) on the right-hand side of Eq.(39) then can 
be extracted by using the spin-echo sequence: 

Usj = Ussexp{-inIjy)U^^exp{i7iIjy) 

n 

= expH4^(2aj7,g{l+ E' {2alh,){2atli,) + ...}] (41) 

l=k<l 

where the sums E' with prime symbol run over all indexes except the index 
j. There are only even-body interactions in addition to the unity operator in 
the bracket {} on the right-hand side of Eq.(41). In order to cancel further all 
these even-body interactions but leave only the unity operator in the bracket 
{} one may first make a simple unitary transformation on the unitary oper- 
ator Usj to convert all these even-body interactions into multiple-quantum 
coherence operators and then uses the phase cycling techniques [23, 27] to 
cancel them (see next sections), 

exp{-ieHjQ) = exp{-i^Fjy)Usj ex.p{i^Fjy) 

n 

= exp(-i4^(2a|7,,){l+ E' (247fe,)(2af7z,) 

l=k<l 

n 

+ E' {2allk,){2atli,){2a;ip,){2allg,) + ...}) (42) 

l=k<l<p<q 

where the eff'cctivc Hamiltonian Hjq is given by 
^77,Q = exp(-^fF,g{^^77,, 

- exp{-i7rIjy)6Hss exp{i7rljy)} exp{i^Fjy) (43) 

and it satisfies 

0HjQ = exp{-qFjy){NeE,, 

— exp{—i'}iFy)N9Ess exp^inFy)} exp{i^Fjy) 

— exp{-i^Fjy) exp{-i7iIjy){N9Ess 

— exp{—i'jTFy)N9Ess exp{inFy)} exp{iT:Ijy) exp{i^Fjy) 

n 

^ W{2a'jlj,){l+ E' {2alh,){2atli,) 

l=k<l 

n 

+ E' (2a|7fe,)(2af7,,)(2a;7p,)(2a^7,,) + ...} (44) 

l=k<l<p<q 
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where the operator Fjp — J2 hp {p — x,y,z). Note that there are the 

identities for the selective phase-shift operation CsiO): 

Cs{±Ne) = exp(Ti^2"£;,,) 

= cxp[Ti(2"^mod(27r))^J = C,(±A^e mod(27r)). (45) 
It follows from Eqs.(41)-(44) that the unitary operator exp{—i9HjQ) can be 
expressed explicitly as 

exp{-i6HjQ) = exp{-i6Hjq) 

= ex.p{-qFjy)Cs{N9 mod(27r)) exp(-i7rFjC,(-A^^ mod(27r)) 
xexp{-i7rIjy)CsiN9 mod{2n)) eicp{inFy)Cs{-N9 mod(27r)) 
X exp{inljy) exp{i^Fjy) (46) 
where the effective Hamiltonian Hj^ is written as 

9Hjq = exp{-i^Fjy){N9 mod{27T)Ess 

— ex-p{—i'KFy)N9 mod(27r)£^ss ex-p{i'KFy)} exp{i^Fjy) 

— exp{—i^Fjy) exp{—i7iIjy){N9 mod{2ir) Ess 

— exp{—i7rFy)N9 mod{27i)Ess cxp{i7TFy)} exp^inljy) exp{i^Fjy). (47) 
One can see from Eqs.(44) and (47) that the effective Hamiltonian 9HjQ 
is proportional to N9 but the Hamiltonian 9Hjq is bounded by the norm 
||^-f^i<j|| < Stt for any qubit number n. It is well known that the highest order 
of multiple-quantum coherence is n for a coupled spin (1=1/2) system with 
n two-state particles [23]. This indicates that the operator Hjq of Eq.(44) 
consists of a variety of multiple-quantum coherence operators whose quan- 
tum orders take any values from — (n — 1) to (n — 1) due to the fact that 
the operator {2ajljz) in the operator Hjq is the LOMSO operator, which is 
also a zero-quantum operator [18]. To cancel all these multiple-quantum co- 
herence operators but leave only the desired term {2ajljz) on the right-hand 
side of the operator Hjq of Eq.(44) the discrete Fourier analysis [28] and the 
phase cycling technique [23, 27] will be exploited below. 

7.1 The discrete Fourier analysis and the phase cycUng technique 

The phase cycling technique is one of the most useful methods in nuclear 
magnetic resonance (NMR) spectroscopy [23, 27]. It has been used exten- 
sively to select the specific multiple-quantum coherences with the desired 
quantum order and cancel any other undesired multiple-quantum coherences 
in the NMR experiments. In principle, the phase cycling technique is based 
on the discrete Fourier transform [28]. The principle of the phase cycling 
technique can be outlined below. Consider a coupled spin (1=1/2) system 
with n nuclear spins. The density operator p{t) of the system then can be 
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classified generally according to different quantum order p 

Pit) = E m- (48) 

p=—n 

A p-order quantum coherence has an important characteristic transformation 
behavior under the z-axis rotations [23, 27] 

exp{—iLpFz)p^{t)exp{iLpFz) = p''(t) cxp(— ipy?). (49) 
That is, a p-order quantum coherence generates a phase shift proportional 
to its quantum order p under the z-axis rotations. This key property is the 
base to separate different order quantum coherences, to cancel the unde- 
sired order and to select the desired order quantum coherences in a coupled 
multi-spin (1=1/2) system. Now by making a series of the z-axis rotations 
with systematic increments of the phase angle ipf. on the p-order quantum 
coherence, 

(Pk^k2n/N, A; = 0, 1,2, ...,Ar- 1, 
and then summing up all the rotational results one obtains 

N-l 

, Af-1 

= />^(t){^ E exp{-i27rpk/N)}. (50) 

fe=0 

Note that there is the exponential sum relation in the conventional discrete 
Fourier transform [28], 

^ e' exp(-i27rp/c/7V) = { J' ^ = ° (51) 

fe=0 P T~ ^ 

where p is an integer and N > \p\. One can reduce the sum of Eq.(50) to the 
form 

^ e' eM-iVkF.)mew{iVkF.) = { (52) 

k=0 P 7^^ 

Now applying the phase cycling technique to the density operator of Eq.(48) 
one obtains with the help of Eq.(52) 

N-l n 

^ E exp(-i(^,F,) E pP(t)exp(z(^,F,) = pO(i) (53) 

k=0 p=—n 

where N > n. Equation (53) shows that by A'"-step phase cycling systemat- 
ically only the zero-quantum coherence components {p = 0) in the density 
operator of Eq.(48) keep unchanged but all other nonzero order quantum 
coherences (p 7^ 0) are cancelled. 

Because the operator Hjq of Eq.(44) consists of multiple-quantum co- 
herences with different quantum orders, by replacing the density operator 
in Eq.(53) with the operator Hjq and with the help of the phase cycling 
technique all the p-order quantum coherences (p 7^ 0) are cancelled in the 
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operator Hjq and one will obtain only the zero-quantum coherences 

N-l 

E eM-i'PjkFjz)HjQeM'^'PjkFjz) ^ H^Q (54) 

fe=0 

where N > n — 1. Note that the operators I^x — ^{Ik + ^k) ^ky = 
^{Ik ~^k)- turns out easily from Eq.(44) and (54) that the zero-quantum 
operator Hjq can be expressed explicitly as 

n 

k^l 

+ E' J«,p,(4+/,+/-/-) + ...} (55) 

where the coefficients Jk,h Jki,pq, etc., depend on the unity-number vector 
{a'l}. The Hermiticity of the operator Hjq shows that the coefficients satisfy 
Jk,i = Ji,kJ Jki,pq = Jpq,ki-i ^^c. Equatiou (54) shows that by iV = n step 
phase cycling one can cancel all the nonzero-order multiple-quantum coher- 
ences in the operator HjQ. However, equations (54) and (55) also show that 
the zero-quantum coherences in the operator HjQ are not cancelled by the 
phase cycling technique. Therefore, one needs to cancel further the residual 
zero-quantum coherences in the zero-quantum operator H^q so as to obtain 
the desired longitudinl magnetizaton term (2a|/j^). 

7.2 The cancellation for the zero-quantum coherences 

There is also the important characteristic transformation behavior under 
the z-axis rotations applied to each individual spin in a spin system, like 
Eq.(49) [23], 

exp{-i2'KpfkIkz)Ik ewii'^T^Pfkhz) = Ik exp(Ti27rp/fc) , (56) 
where fk is called the offset frequency of the kth spin in the spin system. Note 
that one can manipulate each individual spin in the system by an external 
field. By making a series of the z-axis rotations applied to each individ- 
ual spin on the zero-quantum operator Hjq and then summing up all the 
rotational results one obtains from Eq.(55) with the help of Eq.(56) 

iV-l n n 

^ E expH(27rm/iV) E fkIkz]H^Qexp[i{27im/N) E fkhz] ^ 

m=0 k=l,k^j k=l,k^j 

4(2a^^7,-,){l+ E' {JkMlDl^ e' expH(27rm/7V)(/, - /,)]} 



k^l m=0 



n 



+ E {Jkl,pq{h h Ip Iq ) 
k^l^p^q 
N-l 

EexpH(27rm/7V)(/, + /,-/^- /J]} + ...}. (57) 



m=0 
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If all the frequencies of the zero-quantum coherences such as fk — fi, fk + fi~ 
fp — fq, etc. are taken as nonzero integer, it follows from the exponential sum 
relation (51) that all the zero-quantum coherences on the right-hand side of 
Eq.(57) are cancelled exactly, leaving only the desired term 4:{2ajljz) when 
the step number N is larger than the maximum zero-quantum frequency in 
Eq.(57). How to choose explicitly the proper offset frequency set {/i, /2, 
.., /„} so that all the zero-quantum coherences in Eq.(57) can be cancelled? 
The suitable integer set of {/i, /2, .., /„} should satisfy that (a) all the zero- 
quantum frequencies such as fk — fi, fk + fi — fp — fq, etc., take nonzero 
integers and (b) there is a minimum step number N < poly{n) that does 
not divide each zero-quantum frequency. The condition N < poly{n) ensures 
that all the zero-quantum coherences in the operator Hjq can be cancelled 
by the phase cycling of Eq. (57) with a polynomial step number N. It is easy 
to find an offset frequency set {fk} that satisfies only the condition (a). As 
an example, all the zero-quantum frequencies take nonzero integers if the 
offset frequency set {fk} is taken as a super-ascend integer series: fk > 

and fk+i >E ft, e.g., {fk} = {2°, 2^ 2"-^}. However, it may be really 
1=1 

difficult to find the minimum step number A'^ so that N < poly{n) when the 
condition (a) is met. Without losing generality, assume below that all the 
zero-quantum frequencies fzq take nonzero values. With the aid of Eq.(56) 
the following integration identity is constructed similar to the summation of 
Eq.(57) 

T n n 

^ J dt ex.-p[-i27:t E /fe/te]-ff°Q exp[i27ri E fkhz]^ 

-T k=l,kytj k=l,k^j 

n T 

4(2aj/,,){l+E'{4/4+/r)2^ / dtexp[-t{2nt{fk-fi)]} 

k^l -T 

+ E' {JkUni^i-i;) 

k^l^p^q 

/ dtexpH27rt(/fe+/^-^-/,)]}+...}. (58) 

-T 

Note that for a sufficiently large positive number T and any nonzero zero- 
quantum frequency f^q, 

^ / dt eM-i'^^tfzq) = '^^f^r^ = 0. (59) 
-T 27ri fzg 

Inserting Eq.(59) into Eq.(58) and making a variable substitution t — XT one 
obtains 
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, 1 n 

4(2a|/,g = i/rfA{cxpH27rAr E fkhz] 

k=l,k^j 

n 

xH%exp[t27rXT E fkhz] 

k=l,kj^j 

n n 

+ exp[227rAT E fkIkz]H^Qexp[-i2TrXT E /fc4j}. (60) 

k=l,ky^j k=l,ky^j 

The oracle unitary operation exp{—i6ajljz) can be conveniently expressed 
as, by inserting the zero-quantum operator Hjq of Eq.(54) into Eq.(60), 

n-l 

exp(-i^aj/j^) = exp{-«Y6^ E Gxp{-i<^-^,Fj,){HjQ{T) 

+if,Q(-T))exp(?v^^.,,F,,)} (61) 
where the Hermitian operator HjqiT) is written as 

1 n 

Hjq{T) =J dX{exp{-i27rXT E fkhz) 

fe=l,fe7^J 

n 

xi/,Qexp(i27rAr E fkhz)} (62) 

with the matrix element of any pair of the conventional computational base 
|r) and \t) : 

1 n 

{r\H,Q{T)\t) dX{r\H,Q\t) exp[-iXTn E (4 - 4)/^]. 
Obviously, (r|iJjQ(T)|t) = {r\HjQ\t) if the multiple-quantum transition fre- 

n 

quency equals zero, that is, E {'^k~'A^fk = 0; otherwise, {r\HjQ{T)\t) = 

k=l,kj^j 

n 

0. As assumed previously, the zero-quantum transition frequency E 

k=l,kj^j 

{al — a\)fk 7^ where the transition of the pair of the computational base 
|r) and \t) is a zero-quantum transition. Then, {r\HjQ{T)\t) — for any 
zero-quantum transition betwwen the base |r) and \t). For convenient dis- 
cussion, first of all, assume that all multiple-quantum transition frequencies 
including the zero-quantum transition frequencies take nonzero values. Then 

n 

the transition frequencies T[ E (^fe ^ (^k)fk] ^^^^ large numbers, that is, 

k = l,ky^j 

\Tfk\ ^ 1, for a sufficiently large number T, indicating that the integral 
of Eq.(62) contains rapidly oscillating periodic integrand. The conventional 
numerical integration method may not be available for the type of rapidly 
oscillating integrals [29, 30]. This one- dimensional rapidly oscillating integral 
may be first converted into a multiple integral before integrating it numeri- 
cally and this usually will generate an error of 0{j;) [30]. It turns out easily 
that this error is actually zero for the rapidly oscillating periodic integrand 
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operator of Eq.(62) when T is a sufficiently large number. Let = XTf^, 
k = 1,2, ...,n. The one-dimensional rapidly oscillating periodic integral (62) 
then is converted exactly into the multiple integral: 

1 1 n 

HjQ{T) =/ ... / dyidy2...dyn{e-xj£>{-i2'K E Vkhz) 

k=\,k^j 

n 

xHjQ exp{i27r J2 Vkhz)}- (63) 

k=i,k^j 

In fact, the equal matrix element [ifjQ(T)]fc; of any pair of the usual compu- 
tational bases \k) and |/) can be obtained from the multiple integral (63) and 
from the one-dimensional integral (62), respectively. It easily proves from 
Eqs.(44) and (62) that the Hermitian operator HjQ{T) is a diagonal opera- 
tor with the norm ||ifjQ(T)|| = 4 and commutes with the diagonal unitary 
operators exp{±iipji^Fjz). Therefore, equation (61) is reduced to the simple 
form 

eM-iOci^jlj^) = eM-ire^HjQ{T)}eM-ije^HjQ{-T)}. (64) 
This result shows that in the case that all multiple-quantum (including zero- 
quantum) transition frequencies take nonzero values the phase cycling in 
section 7.1 to cancel all the nonzero-order multiple quantum coherences be- 
comes not necessary. There are a number of numerical integration methods 
to calculate a multiple integral [30-36]. One simple numerical method is 
Hua-Wang method [31, 32] based on the number theory [37]. The lattice 
point for the numerical multiple integration is chosen as a real algebraic ir- 
rational point {yi,y2, ...,yn) = u> = (wi, a;2, a;„), where {1,uji^uj2, ■■■,i^n} 
are linearly independent real algebraic numbers over the rationals. Then the 
multiple integral (63) can be replaced with a discrete summation through 

the numerical multidimensional integration [31, 32] 
1 1 

/ ... / dyidy2...dynGjQ{yi,y2, ...,yn) 



1 Mxl 

= 7^n7-l^I ^ ^M,l,k)GjQ{kuj,,kuj2,...,kuJn)+0{M,l) (65) 

[2M + iy k=-Mxi 

where the integrand operator is given by 
GjQ{yi,y2, ■■■,yn) 

n n 

= exp(-i27r E ykIkz)HjQexp{i27r E Vkhz) (66) 

k=l,k^j k=l,k^j 

with the matrix element of any pair of the computational base |r) and \t) : 
{r\GjQ{yi,y2,...,yn)\t) ^ {r\HjQ\t)exp[-in E (4-4)yjt]. 

k=l,k^j 

The integer weight distributions ^{M,l,k) in Eq.(65) are determined from 
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the polynomial identity 

M Mxl 

( E z^Y= E ^M,l,k)z\ (67) 

k=-M k=-Mxl 

and the error function operator 0{M,l) can be derived as [31, 32] (also see 
Appendix B) 

{r\0{M,l)\t) = { ^'^l''^«l'^^(2M+l)sin[7r(m(r,i),a;)]^' ^ (68) 
0, r = t 

where the nonzero integer vector m(r, t) = {mi{r,t),m2{r,t), ...,mn{r,t)) 
with mfe(r, t) = (a^ — a^)/2 = —1, 0, +1; A; = 1, 2, n. The inner product be- 

n 

tween the two vectors m(r, i) and w is defined as (m(r, i), u?) = E <^fe"^fe(^, ^)- 

fe=i 

By using the Chebyshev polynomial of the second kind C/„(cos</?) = sin((n + 
1)(/?)/ sin((/?) the error function of Eq.(68) is reduced to the form 

(r|O(M,0|t) = -(r|//,d0[^|^^ r^t (69) 



where 



;72m(cosv:') 



(2M 

(2M +1) - ^ '^^^'^^ ^ 7r((m(r,i),a>)), {x) denotes the dis- 
tance from the real number x to the nearest integer, i.e, (x) — min({a:} , 1 — 

{x}), and {x} is the decimal part of the real x. By Eq.(44) one can prove 
that \{r\HjQ\t)\ = 4 for any pair of the computational base |r) and |t). Then 
it follows from Eqs.(68) and (69) that the error function operator 0{M, I) of 
Eq.(65) is bounded by 

2" — 12" — 1 

||0(M,/)||<[E E \{r\0{M,l)\t)\y 



r=0 t=Q 

<7^7T^['"e \U,M{coscp{rM'r^'- (70) 

This error function can be made as small as desired when the real algebraic 
irrational lattice point uj and the numbers M and I in the numerical inte- 
gration (65) are chosen suitably. The explicit estimate for the error upper 
bound may be achieved with the Hua-Wang numerical method [31, 32]. As 
shown in Refs. [31-35, 38], for every real algebraic irrational point uj and a 
finite positive number a > 1 there exists a positive constant h — b{a, a;) > 
dependent only on the vector uj and the number a such that 

((m(r, t), u:)) > b{a, u})h{ni{r, t)y (71) 
holds for all nonzero integer vectors m(r, t) (the Schmidt theorem [38]). In the 
Schmidt theorem (71), the distance of the integer vector m(r, t) from the ori- 

n 

gin is defined by /i(m(r, t)) = Yl max(l, \mk{r,t)\). Therefore, /i(m(r, t)) — 1 

k=l 
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for any computational base |r) and \t) in the case of tlie integrand opera- 
tor of Eq.(66). Note that |sin[7r(m(r, t), a;)]| > 2 ((m(r, t), cj)) [31-35]. From 
the inequahties (70) and (71) it follows that the error function operator is 
bounded by 

||0(M, /) II < 4(4^-2'^) V2(^ 3_ y_ (72) 

By choosing suitably the numbers M and /, for example, / ~ n, 2M + 1 ~ 
6(a, w)^^, one can get the desired small error function operator 0{M,l) in 
the numerical integration (65). It is believed that the coefficient 6(a, w) > 
does not decrease exponentially as the dimension size n of the multiple in- 
tegral (63) since it is dependent only on both the algebraic irrational lc 
and the number a [31-35, 38]. As a consequence, the lattice point number 
(2M X / + 1) in the numerical integration (65) does not grow exponentially 
as the dimension size n. 

The Schmidt theorem (71) is still too strong for the numerical multiple 
integration (65), although the integration (65) requires that the Schmidt 
theorem (71) hold only for all those nonzero integer vectors m(r, t) that 
satisfy 

m(r,t) ^ 0, /i(m(r,t)) = 1 (73) 
instead of for all the nonzero integer vectors m(r, t). In fact, the error upper 
bounds (70) and (72) are derived under the previous assumption that all the 
multiple-quantum transition frequencies take nonzero values. This assump- 
tion really corresponds to the condition (73). One can weaken the use of the 
Schmidt theorem (71) in the evaluation of the error function (70) with the 
help of the phase cycling technique in section 7.1. Provided that the phase 
cycling is available in Eq.(61), the Schmidt theorem (71) is required to hold 
only for those nonzero integer vectors m(r, t) that satisfy 

n 

m(r,t) 7^ 0, /i(m(r,t)) = 1, E mk{r,t) - 0. (74) 

k=l 

Obviously, the assumption that all the zero-quantum transition frequen- 
cies take nonzero values corresponds to the condition (74). Although those 
nonzero integer vectors m(r, t) that do not satisfy the condition (74) can 
cause a large error on the right-hand side of Eq.(65), the error is of the 
nonzero-order multiple-quantum coherences and hence can be further can- 
celled by the phase cycling technique, as shown in section 7.1. Thus, inserting 
the numerical integration of Eq.(65) into Eq.(61) one obtains 

n—l 

exp{-iea'jlj^) = exp{-i^ E exp{-i^jk'Fjz) 

k'=0 
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Mxi $(M / A;) 
k=-Mxi [2M + ij' 

+GjQ{-kLJi, -kLJ2, -kuJn)] (iw{'i¥'jk'^jz)-iOo{M, I)} (75) 
where the error function operator is given by 

Oo(M,/) = £ E exp(-i(^^.,,F,,)O(M,0exp(i(^^.fc,F,-,) (76) 

k'=0 

with the error function operator 0{M, I) given by Eq.(68) or (69). Obviously, 
{r\Oo{M,l)\t) = for all the nonzero integer vectors m{r,t) that do not 
satisfy the condition (74). Therefore, one has 

(rlO (M l)\t) = I f conditions (74) hold 

' /I / 1 otherwise ^ ^ 

and the error upper bound is derived as 

l|Qo(M,0||< [ E' |f/2M(cos^(r,t))|^']V2 (78^) 

[2M + ly r,t=0,r^t 

where the sums run only over all pair of the computational base |r) and |t) 
of the zero-quantum transition. It follows from (68)-(71) and Eq.(77) that 
the error function operator Oo(M, /) has the explicit upper bound: 

IIOo(M,oil <W\ii f 

because number of all the linearly independent zero-quantum coherences in 
the quantum system with n two-state particles is given by [23] 

Zo = ( ) - 2". 

Obviously, the error upper bound ||Oo(-^, Oil much smaller than ||0(M, /)||. 

How to express the exponential operator on the right-hand side of Eq.(64) 
or Eq.(75) as a sequence of the nonselective unitary operations and the selec- 
tive phase-shift unitary operations? Choose suitably the algebraic irrational 
point io and the numbers M and I so that the error Oq{M^ I) in Eq.(75) can 
be neglected. Then the Trotter-Suzuki theory [39-41] is exploited to decom- 
pose the exponential operator of Eq.(75) into a sequence of the these unitary 
operations. 

For simplicity, the unitary operation exp(— i^a^/j^) of Eq.(75) is rewritten 

as 

n-l Mxl 

exp{-i9aUjz) ^ exp{-it E E ^jik,k',u})} 

k'=0 k=-Mxl 
n-l Mxl 

xexp{-ii E E Aj{k,k',-u:)} + Oo{M,l) (79) 

k'=0 k=-Mxl 

where t — 9/2 due to the fact that 2ajljz = 1, and the Hermitian operator 
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Aj{k, k', uf) is written as 

xGjQ{kuJi, ku2, kun) exp{iLpjf,,Fj^). (80) 
With the help of the generahzed Trotter formula [39] the unitary operation 
exp{—i9ajljz) of Eq.(79) can be decomposed approximately as 

exp[-ii Aj{k,k',uj)/Lo]}^'> 



exp(—i6a^jL 



n—l 



n-1 Mxl 

{ n n 

k'=0 k=-Mxl 

Mxl 

x{ n n exp[-itAj{k,k',-uj)/Lo]}''° + Oo{M,l)+0{Lo) 

k'=0 k=-Mxl 

where the error function operator 0{Lq) has the upper bound [11]: 

Lo[cxp{—i9ajIjz/ Lq) - 1 + iOajIjz/Lo 
and by Eq.(80) the unitary operator exjp[—itpAj[k, k', u) 
stant p can be explicitly expressed as 
exp[—itpAj{k, k',u:)] = exp{—iLpji.,Fjz) 

, . $(M,Z,A;) ^ , , ,. ^ 

X exp[-2tp ^^^^^^ ^ GjQ{ku;i, k(jj2, fcu;„)J exp{upj^,F, 



(81) 



with any real con- 



(82) 



8n(2M + 1)' - — - -"-v^rjfc'^ J^J- 

The number of the unitary operations exp[—itAj{k, k' ,ijj)/Lq\ in Eq.(81) are 
2Lon(2M x / + 1). It follows from Eq.(81) that the oracle unitary operation 
Uoz{0) is written as 

n 

Uoz{0) =n exp[-2^^a^^/,J 



83) 



n n—l Mxl 

= n [{ n n exp[-it A,(A;,A;',a;)/Lo]}^° 

j=l k'=0 k=-Mxl 
n-1 Mxl 

x{ n n exp[-itAj{k,k',-Lj)/Lo]}^^]+nOoiM,l)+nO{Lo) 

k'=0 k=-Mxl 

where the error function operator nO{Lo) is less than the upper bound [11]: 
max { nLQ[exp{—i6aUjz/ Lo) — I + iOaUjz/ Lq] }. 

This error function can be made as small as desired when Lq is taken as a suffi- 
ciently large number, ensuring that the oracle unitary operation Uoz{0) can be 
expressed correctly as a sequence of the unitary operations exp[—itAj{k, k', uj) 
I Lo] as Eq.(83) within the desired small error (here the error nOo{M, I) can be 
neglected, see Eq.(78b)). Note that Uoy{9) = exp{i^F^)Uoz{0) exp{-i^F^). 
By inserting the oracle unitary operation Uoz{d) {0 = vr/4) of Eq.(83) into 
Eq.(34) one obtains finally the real quantum search network Us (34) which 
is expressed as a sequence of the nonselective unitary operations and the 
selective phase-shift unitary operations Cs{0). 
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In fact, the decomposition (81) for the unitary operation exp{—i9ajljz) is 
the first-order approximated decomposition. A more accurate decomposition, 
that is, the higher order approximated decomposition, may be achieved with 
the help of the Suzuki theory [39-41]. As suggested by Suzuki [39-41], the 
exponential operators on the right-hand side of Eq.(79) is first decomposed 
approximately into a symmetric sequence of simple exponential operators 
exp[—itpAj{k, k',u})] {\p\ < 1) with the error of order 

n-l Mxl 

exp{-it E E Mk,k',<^)} ^ f2m-i({Mk,k',u:)}) + 0(t^-^). (84) 

k'=0 k=-Mxl 

The (2m— l)-order approximated symmetrized decomposition /2m-i in Eq.(84) 
can be constructed as 

f2m-li{Mk,k',U})}) = n S{itp2m-lf) (85) 

j'=l 

with the 2-order symmetrized decomposition S{it): 

S{it) = exp[-itAj{-M x /,0,a;)/2] cxp[~itAj{M x /,n-2,uj)/2] 

X cxp[-itA^{M xl,n- l,Lj)]exp[-itAj{M x /, n - 2, a>)/2] 

X exp[-it A j l- Mx I, Q,uj)/ 2] (86) 

R 

where the parameters {p2m-ij'} are normalized, E P2m-ij' — 1; and it has 

been shown that the 2m-order symmetrized decomposition really equals to 
the (2m — l)-order one [39]. The explicit determination for the parameters 
{P2m-ij'} in Eq.(85) is given in Refs.[39, 40]. By exploiting the generalized 
Trotter-Suzuki formula [39, 40] one can obtain more accurate symmetrized 
decomposition with the smaller error of order 0{L-'-^'^-^h^'^) from Eq.(84) 

n-l Mxl 

exp{-it E E A,{k,k',u;)} = [hm-i{{iA,{k,k',iv)})]^ 

k'=0 k=-Mxl 

+0(t2-/L(2™-i)). (87) 
There are [2n{2M x Z -|- 1) — 1] unitary operators exp[—itpAj{k,k',uj)] for 
each two-order symmetrized decomposition S{it) of Eq.(86). Then there are 
i?[2n(2M X /-I- 1) ~ 1] unitary operators exp[—itpAj{k, k',u>)] in the (2m — 1)- 
order symmetrized decomposition f2m-i{{Aj{k, k' ,uj)}) of Eq.(85) and the 
symmetrized decomposition of Eq.(87) contains the unitary operators with 
the total number of LR[2n{2M xl + l)-l\. 

Now inserting the symmetrized decomposition (87) into Eq.(79) the uni- 
tary operator exp(— i^a^/j^) is expressed as 

exp(-z^^a|/,,) = [/2™_i({iA,(A;, k\ u;)})]^[/2^_i({i^(A;, k\ -a;)})]^ 

+Oo(M,0+0(t2"^/L(2™-i)). (88) 
As the final result, the oracle unitary operation Uoz{0) can be expressed as a 
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sequence of the unitary operations exp[—itpAj{k, k', u?)] which consist of the 
selective phase-shift operatons and the nonselective unitary operations, 

n 

UoM=U {[/2^-l({H(^'^''^)})]^[/2.n-l({H(^,^',-'^)})]'^} 

+nOo(M, I) + nO(i2™/L(2"»-i)). (89) 
It has been shown that the decomposition of Eq.(89) is much more ac- 
curate than that of Eq.(83) and thus, the number of the exponential uni- 
tary operations exp[—itpAj{k,k',u})] in Eq.(89) is much less than that one 
in Eq.(83) with the same error in magnitude [39]. Now number of the 
selective phase-shift operations Cs{0) to compose the oracle unitary oper- 
ation Uoz{d) can be evaluated from the expression (89). The oracle uni- 
tary operation Uoz{0) contains 2nLR[2n{2M x / + 1) — 1] unitary operations 
exp[—itpAj{k,k',u))], whereas equations (46), (66), and (82) show that each 
unitary operation exp[—itpAj{k, k', u))] contains four selective phase-shift op- 
erations Cs{9). Therefore, the oracle unitary operation Uoz{0) really consists 
of the number 8nLR[2n{2M x / -|- 1) — 1] of the selective phase-shift operations 
Cs{6) in addition to the nonselective unitary operations. On the other hand, 
it can be seen from Eqs.(46), (66), (80), (82)-(89) that all the nonselective 
unitary operations including exp{±inFy), exp{±i^Fjy), exp{±i(pjf^Fjz), and 

n 

exp(±i27r J2 Vk-^kz), etc., (note that the index j runs over all n particles 

k = l,ky^j 

in the quantum system so that exp{±i^Fjy), etc., become the nonselective 
unitary operations) in the oracle unitary operation Uoz{0) are one-qubit uni- 
tary operations and the total number of these nonselective unitary operations 
is proportional to the number nLR[2n{2M x Z -|- 1) — 1]. Therefore, the oracle 
unitary operation Uoz{d) can be really expressed explicitly as a sequence of 
the number 0{nLR\2n{2M + — 1]) of the selective phase-shift operations 
Cs{0) and the number 0{nLR^n{2M x / + 1) — 1]) of the nonselective unitary 
operations. Obviously, the quantum network Us of Eq.(34) with the oracle 
unitary operation (89) or (83) could be a polynomial-time quantum search 
network only when the lattice point number (2M x / + 1) of the numerical 
integration (65) docs not increase exponentially as the dimensional size n of 
the multiple integral (63). 

8. Discussion 

As shown in previous sections, the multiple-quantum operator algebra 
formalism has been exploited to construct explicitly a real quantum search 
algorithm. In an unsorted search problem the initial state, i.e., the superpo- 
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sition, and the final state, that is, the target state, are usually given and fixed 
in a quantum system. Then the propagator and its corresponding effective 
Hamiltonian can be constructed explicitly that describe in quantum mechan- 
ics the time evolution of the quantum system from the initial state to the final 
of the search problem. A real quantum search algorithm could be built up by 
starting out such propagator, although the propagator may usually not be 
compatible with the mathematical structure and principle of the search prob- 
lem and hence is not a real quantum search network. There are two families 
of elementary unitary operations, that is, the nonselective unitary operations 
and the oracle unitary operations, for example, the selective phase-shift op- 
erations, in the unsorted quantum search problem. These elementary unitary 
operations are compatible with the mathematical structure and principle of 
the search problem. Then, the propagator is compatible with the mathe- 
matical structure and principle of the search problem and becomes a real 
quantum search network when it is expressed explicitly as a sequence of the 
nonselective and the selective unitary operations. The multiple-quantum op- 
erator algebra formalism plays an important role in the general construction 
of the quantum search algorithm. In particular, the discrete Fourier analysis 
and the phase cycling technique based on the characteristic transformation 
behavior of multiple-quantum coherence operators under the z-axis rotations 
are very helpful for the construction of the quantum search networks. 

An unsorted search problem in a large unsorted database is a hard prob- 
lem in classical computation and there has not been any efficient classical 
search algorithm to solve the NP problem in polynomial time so far. Grover 
has showed that the search problem can be fast solved by a quadratically 
speed-up quantum search algorithm over any classical algorithms [14] . How- 
ever, the Grover algorithm is not really a polynomial-time quantum search 
algorithm. In the paper two possible schemes are proposed to solve the un- 
sorted search problem. One scheme is based on the NMR devices [24] (see 
Appendix A and C). The oracle unitary operations UopiO) {p = x,y,z) in 
the quantum network Us oi Eq.(34) could be implemented directly on the 
NMR devices and hence the quantum network Us becomes a real unsorted 
quantum search network. With the help of the NMR device [24] the unsorted 
search problem which is an NP-problem is converted into a special knapsack 
problem which can be solved efficiently in polynomial time (see Appendix 
A), indicating that the unsorted search problem could be solved efficiently 
on the NMR quantum computer in polynomial time. Another is that the or- 
acle unitary operations Uqp{9) are expressed as a sequence of the nonselective 
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unitary operations and the selective phase-shift operations which can be im- 
plemented directly on an oracle universal quantum computer with the help of 
the phase cycling technique, the numerical multidimensional integration and 
the Trotter-Suzuki theory. It has been shown that the computational com- 
plexity of the unsorted quantum search algorithm is dependent mainly upon 
that of the numerical multidimensional integration. The proper numerical 
multidimensional integration methods should satisfy the requirement that 
the lattice point number to numerically integrate the multidimensional in- 
tegral (63) does not increase exponentially as the dimensional size of the 
multiple integral within the desired error so that the constructed quantum 
search algorithm becomes an efficient algorithm. One of the possible numer- 
ical methods of multiple integration [30-37] to evaluate the multiple integral 
(63) may be the Hua-Wang number-theoretic method [31, 32]. It is believed 
that with the Hua-Wang method the lattice point number of numerical inte- 
gration does not increase exponentially as the dimensional size of the multiple 
integral (63). As a consequence, it is believed that the quantum network Ug 
of Eq.(34) could be an efficient quantum search network on an oracle uni- 
versal quantum computer. Therefore, it is believed that quantum computers 
could solve efficiently a general NP-problem in polynomial time. 
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Appendix A 

An NMR model based on the spectral labeling [24] is proposed to solve 
efficiently the unsorted search problem experimentally. It is based on the 
massive parallelism of quantum computation and the noncoUapse, nonde- 
mohtion, and phase-sensitive measurement in NMR techniques. This NMR 
model could offer the possibility to transform the hard NP-problcm of the 
unsorted search problem into the polynomial-time problem, resulting in that 
the hard problem could be solved efficiently. This result supports the belief 
that a quantum computer could solve in principle a general NP-problem, 
although the present model could work only on a few-qubit system. The 
spectral labeling on the NMR quantum computing is first proposed by Madi, 
Bruschweiler, and Ernst [24] . One of the advantages of the method over the 
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other state labelings is that each resonance peak of the NMR spectrum of the 
ancillary spin corresponds one-to-one to a quantum state of the work space 
of quantum computation. This labeling method requires that the ancillary 
spin be coupled with all the spins in the work space and all single-quantum 
transitions of the ancillary spin be nondegenerate. This may limit the practi- 
cal application of the model in an NMR system with many qubits. However, 
it provides a very simple experiment model to solve efficiently the unsorted 
search problem. 

Consider the weakly coupled spin (1=1/2) system SAMX... as a quantum 
computation device, that is, an NMR quantum computer. The spin S is 
the ancillary qubit and the spins A, M, X, ... form the work space of the 
quantum computation. In particular, assume that all the spins A, M, X, ... 
couple with the spin S. In order to exploit the massive parallelism of quantum 
computation the system is first prepared at a superposition including the 
marked state \s): 

\^) = \r,S) = E'a.|x)[^(|0)-|l))] 

where the ancillary qubit S is at the superposition :^(|0) — |1)). The eval- 
uation of the function f{x) then can be achieved by performing the oracle 
unitary operation Uf on the superposition 

Uf : 1^) - a.k)[^(|Oe/(a;)) - |ie/(a;)))] 

= a.k>[^(|0>-|l»]-a,|.)[^(|0)-|l))] 

x=0,x^s 

where it has been introduced the fact that the function f{s) — 1 and /(r) = 0, 
r ^ s . Obviously, only the target state |s)[:^(|0) — |1))] is inverted and any 
other states keep unchanged when performing each evaluation of the function 
f{x). In general, it is difficult to distinguish the inverted target state from 
the other states in a quantum system after performing only once evaluation 
of the function j{x). This is because each state has the equal probabihty. 
However, there is a significant phase difference between the target state and 
any other state after implementing once evaluation of the function ]{x). 
The phase of the target state is opposite to all of the other states. In the 
Grover quantum search algorithm [14] this phase difference is transferred 
into an amplitude difference between the target state and the other states by 
making the diffusion transformation so that the amplitude of the target state 
is amplified, resulting in that the search is quadratically speeded up. It is 
well-known that the phase difference among states in a nuclear spin ensemble 
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may be detected by the nuclear magnetic resonance (NMR) measurement 
techniques. The NMR signals and spectra carry the information of the phase 
difference. Now, it is possible to reveal the effect of the phase difference on 
the NMR spin ensemble, and especially on the phase and amplitude of the 
NMR spectra if the evaluation of the function f{x) is performed on an NMR 
quantum computer. This phase difference may result in the phase-inversion 
spectrum of the target state with respect to those peaks of the other states. 
As an example. Figure One (see central spectrum) is the conventional NMR 
spectrum of the ancillary spin S* in a four-qubit spin system SAMX. When the 
evaluation of the function f{s) = \a/3/3)) is carried out the target state 
\aPP)\S) is inverted in phase. Now, if the NMR spectrum of the ancillary 
qubit S is recorded one may find the phase- inversion peak of the target 
state (see bottom spectrum). It shows how the resonant peak of the target 
state with inversion phase may be recognized from the others peaks. The 
resonant frequency of the phase- inversion peak can be measured accurately if 
the signal-to-noise ratio of the peak is high enough. The resonant frequency 
actually carries some information of the target state. It is expected to extract 
the information from the resonant frequency. 

The spin Hamiltonian of the weakly coupled spin (1=1/2) system SAMX... 
can be written generally as 

n n n 

H — QsSz+ Yl ^khz + Sz Yl Jskhz+ Yl Juhzhz (^1) 

fc=l fc=l l>k=l 

where the contribution of the decoherence has been neglected; the symbol 
I denotes the spins A, M, X, and Q5, are the chemical shifts of the 

spin S and the A;th spin I, respectively; Jsk is the scalar coupling constant 
between the spin S and the A;th spin I, J^/ the scalar coupling constant of the 
/cth and the /th spin I. Obviously, the conventional quantum computational 
base \r,S) = \r)\S) are the eigenbase of the spin Hamiltonian (Al) with the 
corresponding energy eigenvalue E{r, S), 

H\r)\S)^E{r,S)\r)\S), 

E{r,S)^\bs{ns+ E lalJsk)+ E \alQk+ E Ha^Jki (A2) 

k=l k=l l>k=l 

where the unity-number representation of the eigenbase |r, S) has been used 

(see Eq.(18)), 

|r, S) = (iTi + a[Si) ®{\T2 + a^^S^) ® ... {8)(|T„ + a^Sn) ®{\Ts + hsSs). 
For an arbitrary computational basis |r) of spins / the transition frequency 
of the ancillary spin S is written as 

us{r) = E{r, S = +l/2)-E{r, S = -1/2) = ns+ E H^Sk- {A3) 
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The scalar coupling constants Jsk and the chemical shift are usually fixed 
for a given weakly coupled spin system. If one can measure exactly the 
resonant frequency ^^^(r) in an NMR experiment the unity- number vector 
{a^} will be determined from the above equation (A3) which can be reduced 
to the form 

n 

E hlJsk = fs{r) (A4) 

k=l 

where hi = 2« + l) = 0, 1; A: = 1,2, ..,n and /^(r) = uJs{r)-ns+ E a'^^fc- 

k=l 

Obviously, it is the famous knapsack problem to solve exactly equation (A4). 
It is well known that the knapsack problem is generally an NP-complete 
problem [42, 43]. Then it is usually hard to solve equation (A4). The degree 
of difficulty of the problem (A4) is crucially dependent upon the choice of the 
coefficients of Eq.(A4), i.e., the scalar coupling constants {Jsk}- For exam- 
ple, if the scalar coupling constant set {Jsk, A; = 1, 2, n} is a superascend 
sequence, that is, Jsk > (A; = 1,2, ...,n) and 

k 

Jsk+i >J2 Jsi, {I <k <n - 1), 
1=1 

equation (A4) can be solved efficiently in polynomial time [42, 43]. Once 
the unity number vector {a^} is determined, the oracle unitary operation 

n 

Uop{0) = n exp[— i^a^Jfep] can be directly implemented experimentally on a 

k=l 

universal quantum computer. Then the target state \s) can be obtained by 
making directly the quantum search network Us of Eq.(34) on the superpo- 
sition 1^) in Eq.(19). It must be pointed out that there are a number of 
choices of the coefficients in Eq.(A4) besides the above superascend sequence 
so that equation (A4) can be solved efficiently in polynomial time [42-44]. 

The above NMR model [24] transforms really the NP-problem, i.e., the 
unsorted quantum search problem, to the polynomial-time problem, i.e., the 
special knapsack problem [42, 43]. 

Caption of Figure 1. A simple NMR device with the weakly coupled four- 
spin (1=1/2) system SAMX to solve efficiently the unsorted search problem. 
The spin S acts as the auxiliary qubit. Note that the peak of the target state 
\a(3(3)\S) is inverted in phase with respect to other peaks, as can be seen in 
bottom spectrum. 

Appendix B 
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If the multidimensional function f{xi,X2, ■■■,Xn) can be expanded as the 
absolutely convergent Fourier series: 

oo 

/(xi, X2, = E C(mi,m2, ...m„) 

mi,...mn=—oo 

X exp[—i27r{miXi + m2X2 + ... + mnXn)] (-B1) 
with the complex expansion coefficients C{mi,m2, ■■■rrin) and integer indexes 
{ruk] k = 1,2, ■■■,n}, and if the integer weight functions $(M, are defined 
by the following identity 

M Mxl 

( E z^y= E HM,l,j)Z\ {B2) 

j=-M j=-Mxl 

then one has the numerical multiple integration formula: 
1 1 

/ ... / dXi...dXnf{xi, X2, .., Xn) 


1 Mxl 

= i^M^VM ^ HM,l,j)fiju) + 0{M,l) (53) 
{2M + Ij' j=-Mxi 

with the n-dimensional real algebraic number lattice point uj = {u!i,u!2, ■■,u!n), 
and the error function operator 0{M,l) can be expressed as 

0(M,/) = - S' C(m..m.....m„)( ;fJ^^^+.^),-'"'-")i )- (54) 

mi,...m„=-oo (2M + 1) sm[7r(m, wjj 

where the sums E' with the prime symbol do not include the term with mi = 
1712 = ■■■ = rrin = and the dot product (m, uj) = rriiUJi + 1712^2 + ••• +'^n'-^n- 

Proof: The detailed proof for the above theorem can be seen in references 
[31, 32]. Note that the multiple integral on the left-hand side of Eq.(B3) 

equals the coefficient C(mi, m2, ■■■mn) with mi = m2 = ... = m„ = 0, 

1 1 

C(0, 0, ...0) =/ ... / dXi...dXnf{xi, X2, .., Xn). 


Because the Fourier series of Eq.(Bl) is absolutely convergent one has 

1 Mxl I 00 

f^u^^V ^ HM,l,j)f{jaj)= E C(mi,m2,...m„) 

[2M + ly j=-Mxl [2M + Ij' mi,..m„=-oo 

Mxl 

X E ^{M,l,j)exp[-i2'n-j{miu;i+m2u;2 + — + mnUJn)] 

j=—Mxl 

]^ 00 

^ C(0,0,...,0) + E' C(mi,m2,...m„) 

^^ZiKZ "r J- j mi,...mn=— 00 

M 

X ( E exp[-i27rj("T.ia;i + m2a;2 + ■■■ + rrinOUn)])'' 

j=-M 
1 1 

=J ... J dXi...dXnf{xi, X2, .., Xn) 
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mi,...m„=-oo (zikz + Ij sin[7r(m, cjjj 

where the identity (B2) has been introduced. Obviously, if the error function 
0{M,l) is given by Eq.(B4) then the numerical integration formula (B3) is 
obtained. 

To calculate explicitly the error function operator 0{M, I) in the numeri- 
cal integration (65) the integrand operator (66) of the multiple integral (63) 
is first expanded under any pair of the usual computational bases |r) and \t) 
and the corresponding matrix element then is written as 

n 

{r\Gjg{yi,y2,...,yn)\t) = {r\eHjy\t)cxp[-i2n E mk{r;t)yk] (55) 

where the integer vector m = {mfc(r, t)} with mk{r,t) = — 1,0, +1; k = 
1,2, ...,n for any pair of bases |r) and \t), and mi(r, t) = m2{r,t) — ... — 
mn{r,t) — when |r) = \t). Then, the error function operator 0{M,l) in 
the numerical integration (65) with the integrand (B5) can be found from 
Eq.(B4), 

Obviously, (r|0(M, l)\t) — for any pair of bases |r) = \t). 



Appendix C 

In appendix A the NMR device [24] is used to measure experimentally the 
unity-number vector {a^} in polynomial time. However, this device is quite 
limited in practice. Therefore, one hopes naturally that there is a convenient 
NMR device to measure efficiently the unity-number vector. It had better 
start at the thermal equilibrium state of an NMR quantum ensemble instead 
of the effective pure state [46]. That is, this NMR device can exploit the 
initial mix state of a spin system, for example, the thermal equilibrium state 
to determine experimentally the unity-number vector and hence it may be 
useful in practice. The present device is still based on the massive parallelism 
of quantum computation and the phase-sensitive NMR measurement. 

In section 4 it has been shown that the selective phase-shift operation 
Cs{0) can be equivalent to the oracle unitary operation Uf when the auxiliary 
qubits are used in the implementation of the quantum search problem. For 
example, Cgij) = Uf (equivalent ) when the auxiliary qubit S takes the 

superposition l^*) — — (|0) — |1)) in the superposition: 
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|^i) = |/,5)= E^.|:r)[^(|0)-|l))] (CI) 

x=0 

where the selective phase-shift operation Cs(7r) is apphed only to the work 
qubits /, while the oracle unitary operation Uf, that is, the evaluation of 
function f{s), is applied to all the qubits including both the work qubits 
/ and the auxiliary qubit S. It is also shown that the selective phase-shift 
operation Cs{9) can be expressed as Cs{0) = UfV{9)Uf (equivalent) on the 
superposition: 

N-l 

\^2) = \I,S)= E a^\xm\l) (C2) 

x=0 

where the state of the auxiliary qubits S takes \S) = |0)|1). On the other 
hand, in the matrix representation the superposition (CI) and (C2) can be 
expressed respectively as 

Pi = =Pl7(2)Pl5 (C3) 

N-l 1 

with = ( a^ay\x){y\) and p^g = ^^l^) - + (0|)> and 

P2 = l^2)(^2| =P2/(2)P2S (^4) 

N-l 

with = ( E o-xO-ylx) {y\) and P25 = ® Then there are the 

x,y=0 

following relations when the oracle unitary operations Uf and UfV{6)Uf 
applied to the two superpositions (C3) and (C4), respectively, 

Ufp.uy' = Cs{7i)puC.{7i)-'^p,s, {Cb) 
UfV{e)Ufp,iUfVie)Uf)-' = Cs{e)p,jCs{e)-'^p,s (ce) 

This shows that the action of the selective phase-shift operations Cs{0) and 
Cs(7r) on the work qubits / is equivalent to the action of the oracle unitary 
operations UfV{9)Uf and Uf on the whole system including the auxiliary 
qubits S, respectively. Although the density operators Pi and P2 are pure 
states in Eqs.(C5) and (C6), the two equations (C5) and (C6) still hold even 
when the density operators pu and take any mix states of the work qubits 
/. This is because the unitary operations Uf, V{6), and Cs{6) are linear op- 
erators. Then the equivalent relations: Cs(7r) = Uf and Cs{9) = UfV{9)Uf 
still hold even when the density operators p^j and pg/ a-^® taken as the mix 
states of the work qubits / in Eqs.(C5) and (C6). This point is important 
for the NMR experimental determination for the unity-number vector {a|} 
by startng out the mix state p^^ or of the work qubits /, e.g., the thermal 
equilibrium state, instead of the effective pure state. 

In NMR experiments to determine the unity-number vector {a|} the den- 
sity operators p-^j and P27 ^^-y be prepared as any mix states of the work 
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qubits /, but the auxiliary qubits S must be prepared as the effective pure 
states Pig and p2s, respectively. On the other hand, if the density operators 
Pij and P2J do not include the marked state \s) one has 

Cs{7r)p,jCs{7r)-' = Pu: Cs{e)p2jCs{e)-' = P21. 
In this case it is impossible to determine experimentally the unity-number 
vector {af.}. Thus, it is required that the density operators p^j and P21 be 
taken as any mix state that include any given marked state \s). As a simple 
example, the density operator {pu or P21) including any given marked state 
\s) may be taken as 

n 

pj = aoE+ J2 Skh^i (/i = X, y) (C7) 

fc=i 

This density operator can be generated from the thermal equilibrium state 

n 

of the work qubits: Pq^ = aQE+ ^k^kz by applying a ninety degree pulse. 

k=l 

It is assumed in the following discussion that the auxiliary qubits S are 
always prepared as the effective pure states pig and p2s^ respectively, during 
the action of the oracle unitary operations Uf and UfV{9)Uf so that the two 
oracle unitary operations can be replaced by their corresponding selective 
phase-shift operations Csij^) and Cs{6), respectively, to describe the evolution 
process during the oracle unitary operation 

Pi5 = (^|0)(0| + ^|1)(1| - ^|0)(1| - i|l)(0|) = \e- S^, (C8) 
P2S = |0)(0| |1)(1| = i^E, + Su) (^(^E2 - 

= i^ + i(5l,-^2.)-^1.^2. (C9) 

In order to analyse the evolution process during the oracle unitary operation 
the selective phase-shift operation Cs{0) is decomposed into a sequence of 
elementary propagators built up with the base operators of the LOMSO 
operator subspace [18,19] 

C,(^)=exp(-i^:)expH^:(i: a|27fe,)]expH^:( E alafAhM] 

k=l l>k=l 

n 

xexp[-ie:{ E alatat^SIkJiJn^z)].... (CIO) 

l>k=l 

where 6*g = Og/N. Then it is easy to prove that there are the basic unitary 
transformation relations for the oracle unitary operation 

n n 

Cs{e)h.Cs{e)-^ = h. cos[9:.2{E+ E' apii,+ E' afa^4/;,/^, + ...)] 

(=1 m>l=l 
n n 

+alIkySm[9':.2{E+j:'at2Ii,+ E' afa^47^,7^, + ...)] (Clla) 

1=1 m>l=l 
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n n 

Cs{e)hyCs{e)-' = hyCos[0:.2{E+j:'at2Ii,+ E' afaf^AIiJmz + ...)] 

1=1 m>l=l 
n n 

-aiIk,sm[9:.2{E+j:'at2Ii,+ E' afa^4/^,/^, + ...)] (C116) 

1=1 m>l=l 

where the sums E' I'un over all indexes except the index k. To simplify fur- 
ther the unitary transformation relations (CI la) and (CI lb) the triangular 
functions of the LOMSO operator subspace are expanded in terms of the 
LOMSO base operators [18] 

Ti TL Ti 

cos[9l.2{E+ E' at2Ii,+ E' afa^^iliJmz + •••)] = «o/c^+ E' ^'ki^iz 

1=1 m>l=l 1=1 

n n 

+ E' j'klnP'IlzImz+ E' J'klmp^llzlmzlpz+- (C12a) 

m>l=l p>m>l=l 

n n n ,, 

sin[^:.2(E+ E' af2/^,+ E' afa^4/;,/^, + ...)] = «ofe^+ E' ^kihz 

1=1 m>l=l 1=1 

n I, n „ 

+ E' Jklm'^llzlmz+ E' Jklmp^llzlmzlpz + - (C126) 

m>l=l p>m>l=l 

Now start the initial state pj (C7) {/i = x). By applying the oracle unitary 
operation and then a hard 90° pulse to the spins / the state is transferred 
into 

k=l 1=1 

n ^ n ^ 

+ E Jklm'^^ly^rny+ E Jklmp'^{~^ly^rnylpy) + ■■■} ® P2S 
m>l=l p>m>l=l 
n // " " n II 

+ E £kakIkz{oiQkE+ E' ^ki{-hy)^ T! Jklm'^^ly^my 

k=l 1=1 m>l=l 

~l~ E '^klmp^i. -^ly-^my-^py ) + ...}0P25 (CIS) 
p>rn>l=l 

This state is quite inconvenient for the NMR maesurement of the unity- 
number vector {a|}. However, by applying a purge pulse unit [47] such as 
z-direction gradient field, z-filter, etc., on the density operator (C13) to cancel 
all the multiple-quantum coherences including the zero-quantum coherences 
but leave only the longitudinal magnetization and spin order components 
unchanged, the density operator (C13) is reduced to the form 

Pf = {aQE+ E aQk^kaihz}®P^s {CIA) 
k=i 

The parameters {ckq^.} can be obtained explicitly from Eq.(C12b), 
o^ik = TT E' sin{^^,[l+ E' a|2m,+ E' a|a;4m,mp + ...]} (C15) 

iV mi,mp,mq,... I\ i=l p>i=l 
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where r/ij, nip, rriq,... = 1/2, —1/2 and the sums J2' run over all indexes 
except the index m^. It turns out that the parameters Oq^, = 1/2"^"'^ when 
the phase angle 9s = n/2. Now the density operator pj of Eq.(C14) is quite 
simple. The unity-number vector {a|} can be measured simply from the state 
of Eq.(C14). The measurement may be carried out conveniently by applying 
a hard 90° pulse to the spins / of the density operator of Eq.(C14) and then 
recording the NMR signal of the spins / during decoupling the auxiliary spins 
S. Obviously, the unity-number vector {a|,} may be determined conveniently 
by recording and comparing both two NMR spectra of the density operators 
(C7) and (C14), respectively. 

With respect to the NMR device in Appendix A the present NMR device 
can be even a spin ensemble of linear molecules with neighbor interaction 
and moreover, the auxiliary spins in the device are not required to interact 
with all spins of the work qubits. Therefore, the present NMR device is very 
convenient one to determine the unity-number vector and may be more useful 
in practice. However, the NMR singal of the spins / of the density operator 
(C14) is proportional to l/2"~^, indicating that it decreases exponentially 
as the I-spin number n. Thus, the main drawback of this device is that the 
NMR signal to determine sufficiently the unity-number vector {a^} decreases 
exponentially as the qubit number of the search problem. This is similar to 
the NMR quantum computing based on the effective pure state [48] . 

Although the NMR signal of the density operator (C14) decreases ex- 
ponentially as the qubit number n like the NMR quantum computing on 
the effective pure state, it must be pointed out that they have a significant 
difference. The present device can start at the thermal equilibrium state 
of spin ensemble, while the latter starts at the effective pure state [46, 49]. 
The present scheme to measure the unity-number vector is polynormial-time 
within the realizable size of NMR quantum computation, but the Grover 
algorithm based on the pure quantum state [14] or the effective pure state 
[49] is a quadratically speed-up method even within the realizable size of 
NMR technique. Although in the present NMR device the final NMR signal 
of Eq.(C14) may not be detected due to too low signal-to-noise ratio in a 
spin system with many qubits, this is the NMR technique limit instead of 
the principle limit. Therefore, the final result {a|} will be obtained certainly 
in polynormial time if the sensitivity of the NMR signal of Eq.(C14) is im- 
proved sufficiently in NMR technique [46]. However, in the version of the 
Grover algorithm based on the pure quantum state [14] or the effective pure 
state [49] each quantum state in the superposition has the same probability. 
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then the marked state in the superposition is impossible to be found cer- 
tainly after performing once the oracle unitary operation, i.e., the selective 
phase-shift operation, due to the limit of the quantum measurement principle 
[12] even when the measured signal sensitivity is increased sufficiently. One 
method to find certainly the marked state is to amplify the probability of 
the marked state but suppress all the others. This is just the spirit of the 
Grover algorithm [14]. 

There is an interesting thing. Suppose that the NMR signal of the density 
operator (C14) could be enhanced by applying a sequence of a polynomial 
number of the oracle unitary operations and the nonselective unitary opera- 
tions on the initial state (C7), the reahzable size of the present NMR device 
could be enlarged. 
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